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Abstract 

A combinatorial version of Yamabe flow is presented based on Euclidean triangula- 
tions coming from sphere packings. The evolution of curvature is then derived and 
shown to satisfy a heat equation. The Laplacian in the heat equation is shown to 
be a geometric analogue of the Laplacian of Riemannian geometry, although the 
maximum principle need not hold. It is then shown that if the flow is nonsingular, 
the flow converges to a constant curvature metric. 
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1 Introduction 



In his proof of Andreev's theorem in [29], Thurston introduced a conformal 
geometric structure on two-dimensional simplicial complexes which is an ana- 
logue of a Riemannian metric. He then used a version of curvature to prove 
the existence of circle packings (see also Marden-Rodin [23] for more expo- 
sition). Techniques very similar to elliptic partial differential equation tech- 
niques were used by Y. Colin de Verdiere [7] to study conformal structures 
and circle packings. Cooper and Rivin in [8] then denned a version of scalar 
curvature on three-dimensional simplicial complexes and used it to look at 
rigidity of sphere packings along the lines of Colin de Verdiere. 

Inspired by this work, Chow and Luo [5] denned several combinatorial Ricci 
flows on two-dimensional simplicial complexes, one for each constant curva- 
ture model space. They were able to show that the flows converge to constant 
curvature if a circle packing exists whose nerve is the one-skeleton of the tri- 
angulation. The reader is also directed to some later work of Luo on how 
these flows evolve the conformal structure [20] . We shall use Cooper and 
Rivin's combinatorial scalar curvature to define combinatorial Yamabe flow 
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Fig. 1. Tetrahedron with balls at the vertices. 

on three-dimensional simplicial complexes which is a three-dimensional ana- 
logue of Chow and Luo's work when the triangles are modeled on Euclidean 
triangles. We shall look at the evolution of curvature from a geometric view- 
point, understanding the heat equation on curvature which is induced by the 
flow. The flow turns out not to be parabolic in the usual sense of Laplacians 
on graphs, an analytic property which we study in a related paper [10]. The 
geometric flow perspective is very much inspired by Richard Hamilton's works 
on the Ricci flow, e.g. [14]. 

The combinatorial Yamabe flow is a way of studying prescribed scalar curva- 
ture on simplicial complexes, which we might call the combinatorial Yamabe 
problem. The Yamabe problem has been studied in great detail (see [19] for a 
good overview). The Riemannian case has been solved by the work of Aubin 
[3] and Schoen [28]. Yamabe flow in the smooth category has been studied by 
Hamilton and others. We refer the reader to R. Hamilton [13] and R. Ye [31]. 



2 Geometric structures and combinatorial manifolds 

We essentially take our formalism from Cooper- Rivin in [8]. We shall use the 
notation to denote evaluation of a function / at i in a finite set and / (t) to 
denote evaluation at t in an interval. Let S = {So, Si, . . . , S n } be a simplicial 
complex of dimension n, where Si is the i-dimensional skeleton. We define a 
metric structure as a map 

r : S ^(0, oo) 

such that for every edge € Si between vertices i and j, the length of 

the edge is = r\ + Tj. Any such metric structure, or a particular tetrahe- 
dron within such a structure, is called conformal and the set of all is called 
the conformal class. We can think of this as having an n-dimensional sphere 
packing whose nerve is the collection of edges Si, although it is not necessarily 
an actual sphere packing. One such conformal tetrahedron is shown in Figure 
1. 
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Fig. 2. Tetrahedron with circumscripted sphere. 

Conformal tetrahedra are also called circumscriptible tetrahedra, and the con- 
dition on the edges is equivalent to the condition that there exists a sphere 
tangent to each of the edges of the tetrahedron [1, Chapter 9.B.1] (we call this 
sphere the circumscripted sphere since it is circumscripted by the tetrahedron) 
as seen in Figure 2. 

The function r determines the 2-dimensional faces since there is a one-to-one 
correspondence between triples (n, r^, r^) and triples of sides for Euclidean 
triangles given by 

1 

and so forth. We shall also put the restriction that each higher dimensional 
simplex can be realized as a Euclidean simplex. We shall return to this con- 
dition a little later. 

Each metric of this type is in some sense conformal to the metric where all 
Tx — 1, since they can be gotten by rescaling the function r at each point. This 
is similar to multiplying a Riemannian metric g by a function f 2 at every 
point to get a new metric f 2 g which is conformal to the metric g. The metric 
structure {ri} ieSo determines the geometry, which comes from the lengths of 
the edges, similar to the way a Riemannian metric determines the metric space 
structure of a Riemannian manifold. 

In the sequel we shall limit ourselves primarily to three dimensions. Cooper 
and Rivin [8, Section 3] observe that for a collection {r^, rj, r^, r^} to determine 
a Euclidean tetrahedron, we can use Descartes' circle theorem, also called 
Soddy's theorem, which says that four circles in the plane of radii ri,rj,rk,re 
are externally tangent if 




For a nice proof of Soddy's theorem, see [26]. We also direct the reader to 
the interesting article [18] on Descartes' circle theorem. Looking at the proof 
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Fig. 3. Failure of four circles to be mutually tangent. 



it is clear that if this quantity is negative, then we get three circles which 
are mutually tangent and a circle in the middle which cannot be tangent to 
all the others, as seen in Figure 3, and hence we cannot form a Euclidean 
tetrahedron from spheres of these radii. If Qijke is positive, we can form a 
Euclidean tetrahedron corresponding to {r^, r,-, r k , r e } . So our condition for 
nondegeneracy of the tetrahedron is 




We call Qijke the nondegeneracy quadratic. As noted in [1, Section 793], Qijke 
is really 4 divided by the square of the radius of the sphere tangent to each 
of the edges (the circumscripted sphere), and is related to the volume in the 
following way: 

v ijke = ^ r i r2 j r l r eQijke- 

Thus if we consider the formula for the square of the volume as formal, the 
nondegeneracy condition is that V^ kt > 0. 

Now we shall define a quantity K called the curvature. For a Euclidean tetra- 
hedron with vertices k,£} we define the solid angle a^e at a vertex i as 
the area of the triangle on the unit sphere cut out by the planes determined 
by {i, j, k} , {i, j, £} , {i, k, £} where % is the center of the sphere. Note that the 
solid angle is also sometimes called the trihedral angle. If we define fiijkt as 
the dihedral angle in the tetrahedron {i,j, k,£} along the edge which 
is also an angle of the aforementioned spherical triangle, the Gauss-Bonnet 
theorem gives the formula for the solid angle as 

otijke = Pijke + Pikje + Pujk — 71 '• 

Note that solid angles a^e are symmetric in the last three indices and dihedral 
angles fiijke are symmetric in the first two and in the last two indices. We can 
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now define the curvature Ki at a vertex i as 

Ki = AlX - a ijk£- 

{i,j,k,e}es 3 

Note that the sum is over j, k, £ since the vertex i is fixed. 

The curvature Ki can be thought of as a scalar curvature since it measures 
the difference at a given vertex between the total angles in Euclidean space 
and the total angles of the complex. It was initially looked at by Cooper- 
Rivin in [8]. They found that a metric structure cannot be deformed (staying 
conformal) while keeping the scalar curvature constant. 

Constant scalar curvature is a critical point of the total curvature functional 

T = Y / K t r t . 

Cooper- Rivin showed that the space of nondegenerate simplices is not convex, 
but that the function T is weakly concave as a function of the and strongly 
concave if the condition X r\ = 1 is imposed. We cannot use this to show that 
there is a unique constant scalar curvature metric on a given conformal class, 
but any constant scalar curvature metric is a local minimum of T. 

Another way to prove this is to look at the functional T as a function of 
Si = 1/Vj. Then the set of nondegenerate simplices is convex since it is the 
intersection of a cone with a half-space. However, upon computing the Hes- 
sian of T we can only prove that the function is concave in a neighborhood of 
constant curvature. Much like scalar curvature in the smooth case, there may 
be several constant scalar curvature metrics in a given conformal class. How- 
ever, we have yet to find a complex admitting two constant scalar curvature 
metrics. 

Remark 1 Just because a topological manifold admits a constant sectional 
curvature metric, a given triangulation of that manifold may not admit a 
metric with curvature 0. For any vertex transitive triangulation, i.e. a trian- 
gulation such that the same number of tetrahedra meet at every vertex, the 
triangulation with r\ = 1 for all i e So is constant curvature. If d is the de- 
gree of the vertex, i. e. the number of tetrahedra meeting at each vertex, the 
curvature must be An — d (3 cos" 1 (1/3) — n) . Hence in order for the curvature 
to be zero, d = (Ait) / (3 cos -1 (1/3) — 7r) ~ 22. 795. Thus no vertex transitive 
triangulation of the torus admits a zero curvature metric. This observation is 
similar to the one noted in [2] about dynamical triangulations. 

Remark 2 The curvature considered here is different than the one consid- 
ered in the Regge calculus by Regge [27] and others (for instance, [11], [12], 
[9]). Our curvature is concentrated at the vertices while Regge's curvature is 
concentrated at the edges in dimension 3. The solid angle which we use is 
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reminiscent of the interpretation of Ricci curvature as solid angle deficit (see, 
for instance, the introduction of [4]). 



3 Combinatorial Yamabe flow 



We now define combinatorial Yamabe flow on the metric structure as 

= (!) 
dt y J 

for each % e <So- Note how similar this looks to the Yamabe flow on Riemannian 
manifolds, which is 

where gij is the Riemannian metric and R is its scalar curvature. In particular, 
both preserve their respective conformal class. We use the term 'combinato- 
rial' since this is used by Chow and Luo, but it is really more of a piecewise 
linear or piecewise Euclidean flow because it depends on the geometry of the 
triangulation and not just the topological, or combinatorial, structure. Much 
like the Yamabe flow and Ricci flow, the evolution of curvature will play a key 
role in understanding the behavior of this equation. Next we shall compute 
this evolution. 

Recall the Schlafli formula, which, for a Euclidean tetrahedron denoted by the 
complex {%, T 1: T 2: T 3 } , gives that 

(see Milnor [25] for a proof). We can reorganize this as Cooper- Rivin [8] do 
to get 

n (dp ijke + dpikje + dpie jk ) = J2r-i daiijkt = 

i&% ieT 

or 

da-iju dajike . dakije . daeijk „ 

OTi OTi OTi OTi 

Since there are only four vertices in To, we may denote the solid angle at vertex 
i by «j without fear of confusion. Then we consider 

A=J2 (2) 

ieT 
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so 



dA = E a i d n + E r * den 

— E a i d r i 
ier 

by the Schlafli formula. We thus have 

dA 

dri 

and hence 

dcti d 2 A doij 

dr-j dr-idr-j dr\ 
by commuting the partial derivatives. In our expanded notation, which we 
shall use for complexes larger than one simplex, this says 

da^ke _ dajiu 



dr-j dr\ 

Using our derivation from the Schlafli formula we also have 

dajjke dajjke dajjke dajjke _ „ 
r i k r rj - h rfe - \- re ^ — U. {6) 
dri dr-j dr k dr e 

This equality has a much more geometric interpretation; it says that the di- 
rectional derivative of the angle a^ki in the direction of scaling (rj, rj, r/c, 77) 
is zero, since if these are scaled equally, the new tetrahedron is similar to the 
original and hence all angles remain the same. 

The evolution of curvature is 

d _ d 

j t K * = - E J t <*iM 
az {i,j,k,t}es 3 ai 

y- ( da^ke dr\ | da^edr-j | daijke dr k | da^drA 
{i,j,k7}es 3 \ dr i dr i dr * dr t J 



{i,j,k,e}es 3 



dr-j ' dr k 



+ ^r i (K e -K i ) 
using (3). We call the coefficients 



^_ da^ke 
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In this notation we see that the evolution of curvature is 

d 
di 



-jj_Ki = [^ijke (Kj - Ki) + Q ikje (K k - Ki) + Vt Ujk {Kg - Ki)} . 

{i,j,k,e}es 3 



In order to compute the coefficients fiyH we need to compute the partial 
derivatives of the solid angles. We do this computation using the following 
formulas from Euclidean geometry. Recall that a^kg refers to the solid angle 
of tetrahedron {i,j,k,£} at % and that flij k g refers to the dihedral angle of 
tetrahedron {i,j,k,£} along the edge {i,j}. We also need the face angles. 
Denote the angle of the triangle {i,j, k} at the vertex i by 7^. We can then 
use the law of cosines and the expression for area in terms of sines to compute 
lijk as 

Ijj + ^ik ~ §k _ fj + TjTj + rjfk - r-jr k 
C ° S7ljfc " 2V* " (r t + rj ) {n + r k ) W 

_ 2A ijk _ 2^/r irj r k (n + r 5 + r k ) 
sm y ljk ^ r . + r ^ + r ^ 

where % = r^ + rj is the length of edge {i, j} and Aij k = \Jrirf k {r\ + rj + rg) 
is the area of triangle {i,j, k} by Heron's formula. 

Using the law of cosines for the face angles (4), we can compute the evolution 
of the face angles, which turns out to be 

d Aij k 

where we have introduced the notation P^ k = 2 (r, + rj + r k ) for the perimeter 
of the triangle {i,j, k} . It should be noted that this computation was entirely 
formal, and is thus the same formula derived by Chow and Luo [5] for simplicial 
surfaces. If we define the curvature of a surface to be ki = 2n — X) 7^ then the 
formula for evolution of the face angles implies that the evolution of curvature 
is in fact parabolic in the usual sense of Laplacians on graphs (this is studied 
in greater detail in [10]). The curvature evolution turns out to be 

§ - < A *>. 

where the Laplacian is the one defined by Z. He in [16]. We shall explore this 
aspect in the next section. 

The face angles are used to compute the dihedral angles and solid angles via 
spherical geometry. If we consider the solid angle formed by three planes, say 
those determined by {i, j, k} , {i, j, £} , and {i,k,£}, we see that the planes 
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intersect the sphere and form a spherical triangle. It is clear that the angles of 
this triangle are the dihedral angles flijkt, fiikje, Pujk an d that the length of the 
sides of this triangle are the face angles r yijk,1ije,like, hence the relationship 
between the dihedral angles and the face angles can be expressed in terms of 
the spherical law of cosines, which says 

cos 7^ - cos 7i ife cos -fiji 

coa Pijkt = : — -■ (5) 

sm7 iife sm7y| 



We use formula (5) for the cosine of the dihedral angle and the following 
expression for the volume of simplex {i,j, k,£} 



Vi 



ijki 



2AijkAijt sin Pijkt 



3£{j 



to compute 



dftijke 


2 r . r . r 2 r 2 


1 


l 


dr-i 


'iPijkPijtVijkl 


r l 


rj 








1 




+ (i_r 


)& 


+ — + 




\Tj n, 




rk 



r-j 
2^ 



9f3jjke 
dr k 



r-rfa 



3PijkVij k e [\ri r 



1 

re 

1 1 \ (l 1 1 1 

r k rt. 



Vi \nrk nr e r k r e 



2 + 



Now compute the evolution of the solid angles using the formula for the area 
of a spherical triangle, a ijk t = Pijkt + fiikji + Pitjk - 7r. We get 



da 



ijki 



dr-i 



^ r j r k r t 



^PiikPijlPiktVij 



+ 



Ti \r. 



kt 



1 1 1 
- + — + - 

r k r t 



2 1 1 1 \ 
-+-+-+- 

n r-j r k re) 



+ 



r t, 



rt ( 1 1 1 \ , n 

H 1 1 + (2r; + r-j +r k + r e ) Q ijk e 

n Wi ta r k 



(6) 



which we see is always negative if the tetrahedron is nondegenerate, i.e. Qijke > 
0. The other partial derivatives look like 



da 



ijki 



^ r i r j r k r £ 



dr-j 



-(- + - + -]+-(- + - + - 

; r k re 



"iPijkPijtVijkt \n \rj r k re 
r k re) . 



(7) 



which we would like to say is positive, but is not always (although in the case 
of most "good" tetrahedra, it is positive). 
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Finally we sum cyclically in the last three indices and find 

^ (Pijki + Pikje + Pujk) — ^ijke (Ki — Kj) + Qikji {Ki — K k ) + ^lnjk (Ki — Kg) 



with 



4 2 r 2 2 1 1\ 1/1 1 1 \ 

Q. - H = ^— !LJ — I — I 1 1 -| 1 1 

3P ijk P ije V ijke \n \rj r k r e J rj \r; r k r e J 

( 1 1 



Thus the evolution of curvature is 

dKi ^ / d d d \ 

~dT = ~ { . .fr Vdf ijkl + dt Pikje + Jt^ k ■ 

{i,j,k,t}eS3 \ / 

— E [^ijk£ (Kj — Ki) + Q ikj g (K k — Ki) + Qi£ jk (K e — Ki)] . 
{i,j,k,e}es :i 

We can define our Laplacian as 

(A/) . = ]T (/,- - /,) + - /,) + n tfjfc (/, - /,)] (8) 

{i,j,k,e}es 3 

in order to write the evolution of curvature as 

dK 



dt 



Unfortunately the coefficients flij k £ are not always positive. We notice that A 
is self-adjoint with respect to the inner product 

(f, d) = ]L fi9in 

and satisfies 

E A M = o> 

which is analogous in the smooth category to / Af dx = by the divergence 
theorem. 



4 Combinatorial Laplacians 



In this section we investigate the Laplacian defined in (8) and similar opera- 
tors. Z. He looked at variations of the type we are studying, that is 

dvj 

dt ~ UT ' 1 



10 



Fig. 4. Two triangles with inscribed circles and dual edge. 

in the two dimensions [16]. He then derived a curvature evolution which in- 
volved a combinatorial Laplacian. The work is similar to the smooth derivation 
of the evolution of curvature under a conformal flow on a Riemannian manifold 
given by 

d_ 
d~r 

This was studied by Hamilton in two dimensions for the Ricci flow, and the 
result is 

dR 
~dt 

where A g is the Laplacian with respect to the metric g (see [13]). 



n+dij ~ fdij- 



\f + Rf 



Z. He's Laplacian is the same Laplacian derived by Chow and Luo in [5], which 
can be written as 



using the relation 



A/ 4 = E f(fj-fi) (9) 

*P = v (10) 

where is the angle at vertex % in triangle {i,j,k} and is the length 
of the radius of the circle inscribed in triangle k} . Hence to get (9) we 
simply add 

drj 3 dr-j 3 £ij 

where £*j = r ijk + r ije is the length of the dual edge (see Figure 4). The dual 
vertex k} to the triangle {i,j, k} is the center of the inscribed circle, 

while the edge dual to is the edge which goes from the dual 

vertex {i, j, k} to the dual vertex {i, j, £} which is perpendicular to {i, j} . 
This Laplacian is similar to the Laplacians found in the image processing 
literature which we shall now describe. 

Combinatorial Laplacians on piecewise linear surfaces are used quite a bit in 
image processing, for instance [24], [30]. There is a very clear description by 
Hirani [17] which defines the Laplace-Beltrami operator on functions defined 
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at the vertices as 

1 £*■ 

T/* ^ ' ("^J' ^) 

Vi {i,j}eS! C *3 

where ~jca k is the (n — /c)-dimensional dual of the /c-dimensional simplex a k 
and cr k is the /c-dimensional volume of the simplex. The second line uses our 
notation, where the dual of an edge £*■ and the volume of a vertex V* are 
defined appropriately. Note that Z. He's Laplacian is exactly this, except for 
the volume factor in front, where duality comes from the inscribed circles. We 
also note that in most of the image processing literature, the two-dimensional 
dual comes from the center of the circumscribed circle instead of the inscribed 
circle described here. 

It is interesting to note the geometric justification for the formula (11). If 
we consider the integral of the Laplace-Beltrami operator and use Stokes' 
theorem, we find 

[ AfdV= [ ^-dS 
Ju Jdu an 

where 4^ is the normal derivative and dS is the surface measure. We easily 

an J 

see that if we take U to be it {i} then the normal derivative is 

fj ~ fi 
£■■ ' 

£*j is the surface measure, and V* is the volume measure. 

Recall the definition (8) of the Laplacian we gave in the previous section. This 
Laplacian would be related to the Laplacian in this section (11) if we had an 
analogue of (10). The dual +{i,j} to an edge {i,j} is a surface which goes 
through ^-T for any tetrahedron T containing {i, j} and is perpendicular to 
{i,j} ■ Two pictures of the piece of {i, j} in one tetrahedron can be seen in 
Figure 5. We shall denote the area shown in one tetrahedron by Aijke, where 
the region is perpendicular to the edge {i, j} . The region whose area is AjH 
has four sides: 

• the radius of the circle inscribed in the face k} which intersects 

• the radius of the circle inscribed in the face {i,j, £} which intersects 

• the line from the center of the circumscripted sphere to the face {i,j,k} 
which is perpendicular to the plane determined by {i, j, k} , and 

• the line from the center of the circumscripted sphere to the face {i,j,£} 
which is perpendicular to the plane determined by {i,j,£} . 
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Fig. 5. Two views of the dual to an edge. 

It can also be decomposed as two right triangles, each having hypotenuse 
coming from the radius r^e, of the circumscripted sphere which intersects 
the edge {i,j}- One triangle has a leg coming from the radius of the 
circle inscribed in {i,j,k} and the other has a leg coming from the radius 
rijg of the circle inscribed in {i,j,£}. The legs and hypotenuse meet at the 
same point in {i,j}- We shall call the lengths of the other respective legs 
hijk/ and hijg t k to denote that they are heights of tetrahedra which make 
up the tetrahedron {i,j,k,£}; that is, if we let c denote the center of the 
circumscripted sphere, is the height of tetrahedron {i,j,k,c} with base 
{i,j, k} . Thus the volume of tetrahedron {i,j, k,£} is decomposed as 

Vijkt = X! 7;Aijkhijk,e- (12) 

{i,j,fc}GT 2 6 

If the center of the circumscripted sphere is outside the tetrahedron, we define 
Kjk,e to be negative if the center is on the opposite side of the plane determined 
by {i,j,k} from the tetrahedron. This way, (12) is still satisfied. The hijk,e 
are symmetric in the first three indices, and the fourth indicates to which 
tetrahedron the h corresponds. The area A^t is then computed as 

Aijki — —hijk/Tijk + —hije^Tije, 
where this may be negative. 
We can now show the analogue of (10). 
Lemma 3 

da>ijke Mjkt 
drj £ij 

where is the (signed) area of the dual region to the side {i,j} in the 

tetrahedron {i,j,k,£} as described above. 

PROOF. We simply compute. Recall that is the perimeter of {i,j, k} , 
Aijk is the area of k} , is the radius of the circle inscribed in {i,j, k} , 
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Tiju is the radius of the sphere circumscripted by {i, j, k, £} , and hijk,e is signed 
height of the tetrahedron defined by the center of the circumscripted sphere 
in k, £} and i, j, and k, with base {i, j, k} . The sign of is defined so 
that 

3Vijke = hijk/Aijk + hiji^Aiji + hikejAike + hjke^Ajkt. 
We have the following relations: 

Aijk "^'ijkPijk 

2A{j k ^ j T'iT'jT'k 



Pijk V fi Tj f'k 



We also know that 



4 (I 1 1 1 \ 2 / 1 

— — Qijkt — 1 1 1 ~ 2 ~2 1 2 

rf jke \n Tj r k r e J \rf rj r{ rf ) 

The square of the height hij k / can be computed using the Pythagorean theo- 
rem: 

h 2 - r 2 -r 2 

,l ijk,l ~ ' ijki ' ijk 

= rf jU rf jk fr 1 + r 1 + n _((] L + ]_ + ]_ + ]_Y_ 2 (]_ + ]_ + ]_ + 
4 ^ nrjrk \\ri rj r k r e J \rf r) r\ 

= r ijke r ijk (}_ + J_ + J_ _ J_\ 

4 \ri Tj r k r e J 

In order to determine correctly, we notice that if the center of the cir- 
cumscripted sphere is on face {i, j, k} then 

2 _ 2 
^'ijk ^'ijkl 



and hence 



2 



ji + r. + r,, 1 1 1 IV / 1 1 1 1 \ 
4- 1 k -= (- + - + - + -) -2 +- + - + - 

riffk \n r-j r k r f J \rf rj rf rf j 

(I 1 1 1\ 2 
0= - + - + . 

\n rj r k r e J 

We easily see that if the vector from the center to the side {i,j, k} is in the 
same direction as the outward pointing normal of side {i,j, k} then 

1111 
-+—+—>- 

n rj r k r e 

and if the vector is in the opposite direction than the outward pointing normal 
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then 

1111 
- + — + —<-. 

Ti rj r k r e 

Hence the signed height hij^e is 

, r ijkt r ijk / 1 1 1 1 \ 

n ijk,e — — ^ — — I 1 ) , 

2 \r-i Tj r k r e J 

which may be negative. Notice that this assures that 

3Vijke = hijkjAijk + hij£ t kAiji + h ik £jA ik e + hj k e ti Aj k £, 

where Vij k e is the volume of {i,j, k, £} . The dual area is computed to be 

Aijkt = -hij kt £rij k + -hiji )k riji 

= (\ | 1 | 1 1\ | AA% t fl | 1 1 | i\\ 

4 V p ijk Vi r j r k r e ) Pi je \n r-j r k rj ) 

= gtftu, 1 1 n + L + i \ i n + L + _ (1 _ i ^ 

:5l , ; ;,/ /' ,/,/ ' , / \ r i \ r i r fc r V r i r ^ r ^ ^ r ^ / 

□ 

Corollary 4 We /iawe that 

* {i,j}eSi t% i 

where £*j is the area dual to the side ■ 



PROOF. This follows from the fact that the dual area £*j is simply 

= E A ijk e, 
{i,j,k,e}es 3 

where the sum is over k and £, that is, all tetrahedra incident on the edge 



Note the similarity to Hirani's definition (11). Also note that since £*j may 
be negative, this is not always a Laplacian on graphs in the usual sense (see, 
for instance, [6]). Still, we can prove the maximum principle in more general 
circumstances, as seen in [10]. 
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5 Convergence to constant curvature 



In this section we shall show that if the solution exists for all time, the curva- 
tures converge to a constant. We restrict ourselves to well behaving solutions. 

Definition 5 {r^ (t)} ieSo is a nonsingular solution of the combinatorial Yam- 
abe flow if there exists 5 > such that for each for t e [0, oo) it satisfies (1), 
Qijke > for all k,£} e S 3 , and 

n >5 



for each i G S Q . 

Consider the average scalar curvature 

This can be thought of as an analogue of the average scalar curvature function 

URdv 

f M dV 

on a Riemannian manifold. Note that this curvature really is an average in 
the sense that 

K ■ < h < K 

1 1 mm _ ™ — 1 1 max 

if i^min and K max are the minimal and maximal curvatures. 

Proposition 6 For any nonsingular solution {ri (t)} ieS there is a number 
k (oo) such that the average scalar curvature k (t) and all of the curvatures 
Ki {t) converge to k (oo) as t — > oo. 



PROOF. It is easy to see that k is decreasing along the flow by a direct 
computation: 

dk _ J2i£S Ki r i _|_ (Si€<So Ki r i) 

dt ' T. l&So n (Eies ri) 2 

Rearranging terms we get 

dk _ J2ies ^jes (Ki ~ Kj) r i r j 
dt " (E ie5o r,) 2 
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Furthermore, 



min _ 



where K min is the minimum of the curvature, which is bounded below by 
47r — 2nd ma , x if d ma , x is the maximum number of tetrahedra incident on any 
one vertex. Thus k is decreasing and bounded below, so it must converge to 
a limit k (oo) . Moreover, the time derivative of k must go to zero or k would 
not be bounded below. Hence ^ converges to zero. By formula (13) and the 
fact that 




we see that (Ki — Kj) 2 — > for all pairs of vertices. Hence the curvatures 
becomes constant; this constant must be k (oo) . □ 



6 Long term existence 

In this section we will classify the possible long term behavior. 

Proposition 7 All solutions to the combinatorial Yamabe flow on a maximal 
time interval [0,T) must fit into one of the following categories: 

• It is nonsingular (see Definition 5). 

• T = oo and for some i G S , 



as t — > oo. 

• T < oo and for some k, £} G S3, Qijke — ► as t /* T. 

Solutions with an infinite time interval are covered, and so we need only look 
at finite time singularities. A priori, the following may also happen: 

• there exists i G <S such that r\ (t) — > as t ST 

• there exists i E So such that r\ (t) — > 00 as t S T. 

We can consider Lj (t) = log r\ (t) . Notice that 








dt 



If Li — > ±00 in finite time, then ^j^- 
impossible, though, since 



Ki 



±00 in finite time. This is 



An - 2nd, 



<Ki< 4tt. 
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So finite time singularities occur only because Qiju 0. The case of T = oo 
and there exists % such that 




is somehow analogous to collapse in the smooth case. Hence Proposition 6 is 
an analogue of Hamilton's theorem on nonsingular solutions [15]. 



7 Further Remarks 

It is quite easy to implement the combinatorial Yamabe flow numerically. With 
the help of F.H. Lutz's work (see [22] and [21]) on small triangulations of man- 
ifolds, the author has been able to run examples of the combinatorial Yamabe 
flow on manifolds homeomorphic to the sphere, torus, S 2 x S 1 , and S 2 xS 1 . In 
each case the combinatorial Yamabe flow found a constant curvature metric. 
However, these were all done with relatively small triangulations so many of 
the possible degeneracies coming from a poor triangulation cannot occur. 

As noted in Remark 1, any metric which is vertex transitive, that is, has 
the same number of tetrahedra incident on each vertex, will have a constant 
curvature metric where all are equal. It is not clear that this is the only 
constant curvature metric in the conformal class. Numerical simulation has 
been unable to produce any examples of other constant curvature metrics, 
however. In the case of the torus, the metric where all r\ are equal is usually 
not one where the curvatures are zero. In all simulations of small triangulations 
of the torus tried, the combinatorial Yamabe flow converges to a positively 
curved constant curvature metric. 

Constant curvature as indicated in this paper is a weak condition. A metric 
such that the sum of the dihedral angles along each edge is equal to 2ir is a 
Euclidean structure on a manifold, but a metric where the solid angles equal 
4n is not necessarily. The difference is something like the difference between 
constant sectional curvatures of zero and constant scalar curvature of zero; it 
is possible to have a metric where the latter is true but not the former. 
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